Ferromagnetism of the 1-D Kondo Lattice Model: 



: A Quantum Monte Carlo Study 



3 

Hi 

in 



13 



M. Troyer and D. Wiirtz 
Interdisziplinares Projektzentrum fiir Supercomputing 
and Institut fiir Theoretische Physik, ETH Ziirich, 



> 

o 

i> 

O ; CH-8092 Ziirich, Switzerland 

Abstract 

g 

o 

o 



The one dimensional Kondo lattice model is investigated using Quantum 
Monte Carlo and transfer matrix techniques. In the strong coupling region 
ferromagnetic ordering is found even at large band fillings. In the weak 
coupling region the system shows an RKKY like behavior. 

In recent years the Kondo lattice model (KLM) has been investigated to de- 
scribe the physics of the so called Heavy Fermion systems [|TJ. The question was 
addressed whether this simple model could account for the rich variety of phases 
found in Heavy Fermion materials, paramagnetism, antiferromagnetism as well as 
superconductivity. However, as a typical model of strongly correlated electron sys- 
tems it could be analyzed by only few approximate treatments^, [| and has still 



resisted to give clear insight into the various possible ground state. It would be 
of particular interest to understand the phase diagram of the KLM. In this letter 
we will tackle this problem by using Quantum Monte Carlo techniques for the one 
dimensional KLM. The results found here partially contradict the phase diagrams 
given by Ref. || obtained by variational methods. 

The KLM consists of a lattice of L localized spins (Sij = 1/2; % = 1,...,L) 
coupled to a single band of conduction electrons (creation operator cf CT ; i = 
1, L; a =|, |). It is described by the Hamiltonian 



o-=T,| i=i i=i 

( Ci \ 



2 ( c h c ll ) d 



(2) 



where a are the usual Pauli matrices. The first term denotes the hopping of conduc- 
tion electrons between nearest neighbor sites. The second term is an antiferromag- 
netic (J > 0) exchange coupling between localized spins and conduction electrons. 
The model can be derived from the periodic Anderson model in its strong coupling 
limit §. 

The behavior of the KLM is understood only for some limiting cases. At half 
band filling the ground state is a spin singlet || . The nature of this singlet changes 
from localized singlet pairs in the limit J/t — > oo to a collective singlet at smaller 
values of J / 1. A spin gap has been found at intermediate values of J / 1 using Quan- 
tum Monte Carlo methods [0] and for the entire range of the coupling using exact 
diagonalization || . On the other hand in the case of only one conduction electron 
the ground state of the system is known to be an incompletely saturated ferro- 



magnet |J. The same behavior has been found for two particles and large values 
of J/t J[0|. Questions arising immediately are: what is the ground state of the 
system away from these special cases, where are the boundaries of the ferromag- 
netic region? Contrary to expectations from RKKY mechanisms our results show 
ferromagnetic ordering in the strong coupling region even for very large fillings. 
In the weak coupling limit an RKKY like behavior is observed at intermediate 
temperatures. 

We have employed a generalization of the standard Quantum Monte Carlo 
World Line Algorithm (WLA) and Quantum Transfer Matrix Algorithm (TMA) 



11 . The standard checkerboard WLA has to be extended to add the localized 



spins ||12|| . There is no hopping in the localized band, therefore the only local move 
we have to add is one that can exchange an up (down) spin in the conduction 
band with a down (up) spin in the localized band at any point of the checkerboard. 
This move accounts for the exchange interaction between the two bands. A global 
move that can change the spin of a world line allows for fluctuations in the total 
magnetization. The systematic error due to the finite Trotter time step At is 
controlled by extrapolating from results for Art = 0.25 and Art = 0.5. The usual 
zero winding number boundary conditions are used. 

Due to the fermionic degrees of freedom there is a sign problem. This sign 
problem is most severe at values of J ~ t and near fillings of p = 2/3, where 
p = N/ L and iV is the number of conduction electrons. For larger or smaller values 
of J/t and for a filling of p = 1/3 the sign problem is not too severe to make Monte 
Carlo simulations impossible. 



To study the magnetic properties of the KLM we have measured the charge and 
spin structure factors and susceptibilities. The structure factors are defined as the 
Fourier transforms of the real space correlations. For the conduction band we have 
for the charge and spin structure factors: 

I L 

Scharge{q) = ^ ^ K\t + n j,l) K,T + n Kl) , ( 3 ) 
j,k 

Ss^M) = \ E e^- k) (n 3A - - n jti )(n kA - n Kl ) = ±jh e^S*S z k . (4) 
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In the same manner we define the spin structure factor for the localized spins: 

W9) = xE ei! °' _t) W- ( 5 ) 

j,k 

The static susceptibility can be calculated as 

X (q) = j [drj^e^-V Tr ((5/ + S^e^Sjl + Sye~^ H ) , (6) 

^ - 70 j,k 

where f3 is the inverse temperature 1/T. 

Simulations have been performed on lattices of L = 12, 18 and 24 sites, at band 
fillings of p = 1/3,2/3 and temperatures down to fit = 32. Coupling constants 
J/t — 0.2, 0.5, 1, 2, 4 and 10 have been investigated. Within the error bars of our 
results we cannot see any finite size corrections. From this we conclude that the 
results show properties of the thermodynamic limit. In the figures we only show 
the data for the largest reasonable sizes allowed by the sign problem. 

In the strong coupling region we observe a tendency toward ferromagnetic or- 
dering at both fillings. The spin structure factor of the localized spins and the 
susceptibility both show a clear peak at q = (Fig. 1,2 (a)). Furthermore the 
q = component of the susceptibility is the fastest diverging one when lowering 
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the temperature. We note that this ferromagnetic ordering found at large band 
fillings of p = 1/3 and p = 2/3 is not expected from the RKKY mechanism but is 
a new characteristic of the strong coupling region. Also a slight tendency toward 
ferromagnetic ordering can be observed in the spin structure factor of the conduc- 
tion band (Fig. 2(b)). Double occupancy of a site is strongly suppressed as it 
costs energy of the order of J, leading to an effective on-site repulsion. Therefore 
the charge structure factor is essentially that of spinless fermions, showing a 4fcp 
structure (Fig. 3). 

To get more insight into the magnetic properties we have studied the temper- 
ature dependence of the static uniform susceptibility x{q = 0) f° r J/t — 4 and 
p = 1/3 (Fig. 4). We have simulated systems of L = 6,12,18 and 24 sites by 
the Monte Carlo method and the L = 6 site system with the TMA. Special care 
is necessary for L = 6 because the ground state is completely different depending 
on the boundary condition. It has been shown using exact diagonalization that 
the ground state is a spin singlet for periodic boundary conditions, while it is an 
incompletely saturated ferromagnet for antiperiodic boundary conditions ||10|| . At 
J/t = 4 the ground state energy of the ferromagnetic state is lower than that of the 
singlet state. We have calculated x(T) both for periodic and antiperiodic bound- 
ary conditions in the canonical ensemble using the TMA. In the high temperature 
regime x{T) can be obtained for an infinite size system in the grand canonical en- 
semble by using the M — 1 approximation (M is the Trotter number) in the TMA 

At high temperatures the susceptibility is the sum of that of a free conduction 
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electron system and that of free localized spins Tx = \{1 + p — p 2 /2). When the 
temperature is lowered to about T J the conduction electrons start to lock into 
singlets with localized spins and the susceptibility reduces to about the value for 
the remaining spin degrees of freedom Tx ~ |(1 — p)- Lowering the temperature 
even further the spins start to order ferromagnetically, leading to an increase in 
Tx and possibly a divergence in an infinite size system. 

At small values of J ft a completely different behavior can be observed. The 
effective on-site repulsion is smaller, leading to an increase in the 2k f component 
of the charge structure factor. The charge structure factor resembles that of nearly 
free electrons. At the same time the q = component of the conduction electron 
spin structure factor and susceptibility is reduced while a cusp emerges at 2/cp. In 
the localized spins a 2kp structure is induced by the conduction electrons. This 
behavior may be called RKKY liquid. At lower temperatures the 2kp peaks in the 
spin structure factors become more pronounced while at the same time the q = 
component is suppressed. From this we conclude that in the temperature regime 
of our simulations {j3t ~ 32) the system is dominated by the RKKY interaction. 
It does not show a Heavy Fermion behavior, characterized by a wave vector k' F = 
(p + 1)tt/2, including both conduction and localized electrons. 

A problem arising in the small J / 1 region is that the effective coupling gets very 
small. Although the temperature of the Monte Carlo simulations is well below the 
Fermi temperature, it is still orders of magnitude above the single impurity Kondo 
temperature T K . T K « 1CT 5 (KT 9 ) at J/t = 0.5 and p = 1/3 (2/3). This makes 
it hard to get information on the ground state of the system from Monte Carlo 
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simulations. 

In figure 5 we summarize the results of our Monte Carlo simulations. At large 
values of J/t we find ferromagnetic ordering. There the q = component of 
the susceptibilities is the dominant and fastest diverging one. From our result we 
conclude that the ferromagnetic state is stable for a much wider region of the phase 
diagram than suggested in Ref. ||]. In the small J/t region the system shows an 
RKKY liquid behavior, where the q = 2kp component of the susceptibilities and 
spin structure factors is the dominant one. At J/t = 1 and p = 1/3 the system 



shows a transitional behavior. Probably this point is near the phase boundary 

The calculations were done on the Cray Y-MP/364 of ETH Zurich and on the 
NEC SX-3/22 of the Centro Svizzero di Calcolo Scientifico CSCS Manno. The work 
was supported by the Swiss National Science Foundation under grant number SNF 
21-27894.89 and by an internal grant of ETH Zurich. We want to thank T.M. Rice, 
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Fig. 1 Static uniform susceptibility for the Kondo lattice model at p = 1/3 
and J/t ranging from 0.5 to 4. The temperature is (3t = 24. At this 
filling 2k F = 7r/3. 

Fig. 2 Spin structure factors for (a) the localized spins and (b) the conduction 
band for the same parameters as in Fig. 1. 

Fig. 3 Charge structure factors for the same parameters as in Fig. 1. 

Fig. 4 Temperature dependence of the static uniform susceptibility for a fill- 
ing of p = 1/3. The solid lines show the M = 1 approximation for high 
temperatures and the transfer matrix (TMA) results for periodic (PBC) 
as well as antiperiodic (APBC) boundary conditions. Monte Carlo data 
for L = 6 and L — 12 sites and zero winding boundary conditions are 
included. Due to the sign problem the results for larger lattices have 
much larger error bars. They show no qualitatively different behavior. 

Fig. 5 Phase diagram of the Kondo Lattice Model. The Monte Carlo sim- 
ulations have been performed at the points shown here. The different 
symbols denote points with ferromagnetic, transitional and RKKY liq- 
uid behavior. 



9 



